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The one-dimensional polynuclear growth model with external sources at edges 
is studied. The height fluctuation at the origin is known to be given by either the 

^ . Gaussian, the GUE Tracy- Widom distribution, or certain distributions called GOE^ 

k> , and Fq, depending on the strength of the sources. We generalize these results and 

^ ' show that the scaling limit of the multi-point equal time height fluctuations of the 

model are described by the Fredholm determinant, of which the limiting kernel is 
explicitly obtained. In particular we obtain two new kernels, describing transitions 
between the above one-point distributions. One expresses the transition from the 
GOE^ to the GUE Tracy-Widom distribution or to the Gaussian; the other the 
transition from Fq to the Gaussian. The results specialized to the fluctuation at the 
origin are shown to be equivalent to the previously obtained ones via the Riemann- 
Hilbert method. 
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1 Introduction 

Surface growth is an interesting phenomenon in nature [1]. In particular various shapes 
show up due to the interplay of non-linearity, fluctuation and boundary effects. In [2], 
Kardar, Parisi and Zhang proposed a non-linear stochastic differential equation called 
the KPZ equation. This equation, although now considered to be insufficient for the 
description of most realistic situation in nature, defines a universality class of the surface 
growth phenomenon and plays a prominent role in the theoretical study. 

In the one dimensional case we can analyze the KPZ equation exactly. The rough- 
ness and dynamical exponents were obtained by the renormalization techniques [2] and 
the Bethe Ansatz method [3,4]. Recently we have been obtaining more sophisticated in- 
formation for the height fluctuation in the one-dimensional KPZ universality class; not 
only the exponents but also the height distribution itself have been obtained. For various 
probabilistic models belonging to the KPZ universality class [8-13], it has turned out that 
the height fluctuation is equivalent to the Tracy- Widom distribution [5,6], the distribution 
of the largest eigenvalue in some random matrix ensemble [7] . 

Among these models we focus on the polynuclear growth (PNG) model. For the PNG 
model the relation to the random matrix theory was first pointed out by Prahofer and 
Spohn [14]. They mapped a specific PNG model to the longest increasing subsequence 
problem in random permutations and then applied the Baik-Deift- Johansson theorem [15]. 
The dependence of the height fluctuation on the geometry of the model is also studied 
based on the works of Baik and Rains [16-18] in random permutations with symmetries. 
As a result the deep connection with the various universality classes in random matrix 
theory (RMT) have been revealed. For example, for the droplet growth in an infinite 
line, the height distribution can be described by the GUE Tracy- Widom distribution [14] 
whereas in a half-infinite line, the height fluctuation at the origin can be described by 
the GSE/GOE/Gaussian fluctuation according to the strength of the nucleation rate at 
the origin [19,25]. On the other hand the GOE Tracy- Widom distribution represents the 
height fluctuation in a growth on a flat substrate [19,20]. 

Next we broaden our field of view from the height fluctuation at one point to that over 
some region of the surface. In other words, we address the issue about the roughness of 
the surface. The spatial configuration of the height fluctuation is expected to converge to 
the universal stochastic process after a proper scaling when the space direction is treated 
as the time direction of the process. In general a stochastic process is characterized by a 
dynamical correlation function together with a fluctuation at one point. Thus if we try 
to understand the universal aspect of the spatial configuration of the fluctuation, we need 
information about the multi-point correlations of the height fluctuation. 

The correlation between distinct points of the PNG model is closely related to the 
multi-matrix model. In [21,22], the multi-point equal time height fluctuations of the one- 
dimensional PNG model was studied for a droplet initial condition in an infinite space. 
It was found that the correlation is described by the Airy process, which is the process 
of the largest eigenvalue in the Dyson's Brownian motion model [29] for GUE. This also 
appeared in the facet fluctuation in the crystal [23,24]. In [25], similar quantities are 



evaluated for a special value of the external source in case of droplet growth in half space. 
It was shown that the correlation is described the Dyson's Brownian motion model which 
has transition between GOE/GSE to GUE. The same process also appeared in the vicious 
walk problem [26]. 

In this paper we study the PNG model with external sources. The goal of this study 
is to understand the universal process characterizing the roughness of the surface for this 
model. There are mainly two different regions. One is the region where the edge effects 
are dominant and the other is the region where the bulk dynamics prevails. The statistics 
of height fluctuation obey the one dimensional Brownian motion near edges and the Airy 
process in the bulk. When seen from far away, the above two types of regions are separated 
by a well distinguishable point. But we can focus into a small region around this point, 
in which the edge effects and bulk dynamics are competing. We are especially interested 
in this region because there appear new processes describing the transition of the surface 
fluctuation. 

For this purpose we obtain the Fredholm determinant representation of the multi-point 
correlation function in the bulk region, near edges and around the special points mentioned 
above. We summarize the results as Theorems 3.1, 4.1 and 5.1. In particular the Fredholm 
determinant expressing the correlation in the intermediate region appears in this paper for 
the first time. They describe the transition between GOE^ and GUE/Gaussian, or Fq and 
Gaussian, where GOE^ means the distribution of the larger of the largest eigenvalues of two 
independent GOEs and Fq is a certain probability distribution which has no interpretation 
in RMT [16-18,27]. 

The paper is organized as follows. In the next section, we recall the definition of 
the model and some known facts are reviewed. In section El we give the description of 
the equal-time multi-point fluctuation for finite system. The asymptotic results for fixed 
values of parameters are also stated. In section^ the transition near the GOE^ fluctuation 
is discussed. In the following sectional the transition around Fq is studied. In both section 
13] and El we also discuss the connection between Fredholm representation and the Riemann- 
Hilbert representation of Baik- Rains [27]. The last section is devoted to the conclusion. 

2 Model and One-point Height Fluctuation 

In this article we mainly consider the discrete PNG model with external sources studied 
previously in [19,27]. First of all, we briefly explain the PNG model. The PNG model is 
a simple model of layer- by-layer growth [1]. The discrete version of the model consists of 
the following three rules. 

1. nucleation: A nucleation with height k is generated according to the geometric 
distribution. An object made by this rule is called a step. 

2. lateral growth of a step: Once a step is produced, it grows laterally by one step 
in both directions during each time step. 



3. unification of steps: When two steps with distinct heights collide, the height in 
the colliding point becomes that of the higher step. 

Note that the rule 1 is only probabilistic and the last two rules are deterministic. These 
rules are illustrated in Fig.l. 

We can formulate the above rules of time evolution mathematically as follows. Let r ^% 
and t G N = {0, 1, 2, ■ ■ ■ } denote the discrete space and time coordinates respectively and 
/i(r, t) the height of the surface at position r and at time t. The rules 1 ~ 3 can be collected 
as 

/i(r, t + 1) = max(/i(r - 1, t), /i(r, t), h{r + 1, t)) + cu(r, t + 1), (2.1) 

with the initial condition /i(r, 0) = 0. Here uj is the random variable expressing the height 
of nucleation and takes a value in N. a;(r, t) = if t — r is even or if |r| > t, and 

w{%,2)=uj{%-2,%^2-\), (2.2) 

(i, j) G l\ are geometric random variables. The parameter of this random variable is taken 
to be of the form ajfej, 

^\w{%,3) = A;] = (1 - a,hj){a,hjf, (2.3) 

for A; G N. In order to consider the effect of external sources at both edges, we take 

(2.4) 

(2.5) 

The parameter a is related to the frequency of nucleations in the bulk whereas the pa- 
rameter 7-1- represents the strength of the external sources at the edges. The bigger 7-1- 
is, the stronger the source is. In the following, we assume 7_ > 7-1- when we do some 
computations; the results for the case where 7-|_ > 7_ is obtained from the symmetry. 

Some snapshots of Monte Carlo simulations are given in Fig.2. When the external 
sources are not very strong (Fig. (a)), the effects of the external sources are important for 
some region near the edges whereas the bulk dynamics is important for the curved region at 
bulk. On the other hand, when the external sources are strong (Fig.(b)) they are dominant 
for the whole region. In the critical situation (Fig.(c)), they control the whole region but 
a certain point. More precisely the shape is described as follows. Let us set 

<^^ = Y^ (" + V^^^^) ' (2-6) 

a(l-2a7 + 7^) ^ ^(7^-1) ^ .. 7^ 

"^^(^' ^) = (7-a)(l-a7) (7 - «)(1 - «7)^' ^ ^ ^ 




and 



(l-a^)(7^-l) 
1 + o;^ — 4^7- + 7^ + d^^i 
(l-^^)(7^-l) 
1 + a^ — 4q;7_|_ + 7! + d^^'i 



^ 1 ■ -.,2_/i^,^, 1^,2 , ^,2^,2' (2.^ 



/^+ = - . , J ..,,,. ,0 2 - (2-9) 



Notice /9_ < /3+ -v^ 7+7- < 1- Then the thermodynamic shape is given by the following. 
(i) When (3- < (3+, 

h{r = 2(3N, t = 2N)/N ~ I a(/3), (3. < (3 < f3+, (2.10) 

[aG+(/5,7+), (3>(3+- 

(ii) When (3_> (3+, 

hir = 2pN,t = 2N)/Nr.h-lfl^^-\' ^ < ^- (2.11) 

[aG+(/9,7+), P>Pc, 

with Pc being the solution of aG-(/3,7-) = aG+(/3,7+). 

The fluctuation properties of the model change drastically at the connecting points (3± 
and f3c of the limiting shapes. Let us define the two scaled height variables. The first one 
is 



^ ^ ^^ Aii^r - zpoJV -hzci^PoJJV^r, t zi\ ) aypo ^ ^^/,, )im 


(2.12) 


c^(/3o)iV3 


where 




1 
(1 - d){l - I3^)(> 


(2.13) 


c{l3) = a-5(l - /?2)t(^l + /3 + a^l - /?)5(v/l - /? + a^l + /3)i 


(2.14) 


The second is 





rr(G±),^ ^,_ Hr = 2(3oN, t = 2N)- aG±i/3o, ^)N 

c?g(7)^2 

where 

^^7(1 + a2 _ 4^7 + 72 + a2y) 

^^^^^ = (l-a7)(7-c.) • ^'-''^ 

About the one point height fluctuation, we have the following. 



Theorem 2.1. 

(z) When 13. < (3+. 

a) For 13^ < [3q < (3+, 

lim P[iJ^(0,/3o)<s]=F2(s), (2.17) 

where F2 denotes the GUE Tracy-Widom distribution, which is the distribution of the 
largest eigenvalue in GUE [5]. 

b) For(3o=P- or/3o = (3+, 

lim ¥[Hn{0,(3o) <s]= Fi(s)2. (2.18) 

where Fi is the GOE Tracy-Widom distribution [6]. Thus Ff, which is denoted by GOT?', 
means the distribution of the larger of the largest eigenvalues in two independent GOEs. 

c) For Po < P- or Po > P+, the fluctuation is Gaussian. One has 

lim F[hI^^\Po,1±) <s]= ^ f die-^^^^^. (2.19) 



«? 


«i dG('<-f 


g 2(/3_-/3o) 


r^ p 2(/3o-/3+) dG(7+)2 


v/27r(/3_ - /5o) ^ 


U V2vr(/5o-^+)g^ 



(11) Whenp. > p+. 

a) For Po < Pc or Pq > Pc, the fluctuation is Gaussian (cf. (i-c)). 

b) For Pq = Pc, the fluctuation might be given by 



lim P[i7if-)(/5o,7-) <s]= di,^===== / di, I _ ^\ ^u:r^- (2.20) 

J —00 

(Hi) When /5_ = P+. 

a) For Pq < P- or Pq > P-, the fluctuation is Gaussian (cf. (i-c)). 

b) ForPQ=p., 

lim P[i7^(0,/5o)<s]=Fo(s). (2.21) 

W— >oo 

Here Fq is a certain distribution with mean zero explained in [17, 18,27]. 

These are obtained as corollaries of the theorem in the next section. The special case 
of the height at the origin {Pq = 0) was previously studied in section 4 of [27] using 
the connection of the problem to the combinatorics of Young tableaux. The limiting 
distribution was obtained using the Riemann-Hilbert method and the results were given 
in terms of the solution to the Painleve equation. In [19], for the continuous model, the 
basic picture of this theorem was expected based on a physical argument but has not 
been shown explicitly. These distributions also appear in the fluctuation properties of the 
one-dimensional asymmetric simple exclusion process (ASEP) [9,10]. 



3 Multi-point Height Fluctuation 

As was observed in [21,22,25,28], the equal time multi-point correlation of the PNG model 
can be analyzed by extending the original model to the multi-layer version. The weight of 
the multi-layer version is equivalent to that of non-intersecting many-body random walk. 
In particular, for the model under consideration, this can be borrowed from the results 
in [22] . Following [22] , we only consider an odd time M = 2N — 1 in the sequel. Let us 
consider the weight for the configuration {x^j} = x oi n non-intersecting paths from the 
time r = —M to r = M given by 

M-l 

w„,Af(x) = Yl det (0^,^+1 (x[,a;^+^))"j=i, (3.1) 

r=-M 

where 



^2i-l,2i 



y <x, 



^2,,2,+i(x,i/) = <;;j Jl^^' (3.3) 

L ojv_jjOjy_^-, y ^ X, 

and xf^ = xj^ = 1 — i {i = 1,2, ■ ■■ ,n) is fixed. Note that the same weight gives a 
weight for a time evolution by the time M of the PNG model. For the PNG model with 
external sources we consider in this paper, the parameters aj,b/s are taken to be ()2.4j) . 
()2.5|) . Strictly speaking, the weight of the multi- layer PNG model and the weight ()3.1|) 
with ()2.4|) . (12.51) and n = N are slightly different. Essentially the same remark was already 
given in [25]. The difference is, however, negligible in the scaling limit in which we are 
mainly interested in this paper. 

For each fixed 7±, as iV -^ oo, we have the following results. 

Theorem 3.1. 

(z) When 13 _ < [3+. 

a) For /?_ < /?o < /?+; ^he equal time multi-point distribution function is described by the 

following Fredholm determinant. 



lim W[Hn{Ti,I3o) < Si, ■ ■ ■ , HN{Tm, Po) < Sm] 

TV— >oo 

oo ^ m „ m „ 

= 5Z^ 5Z / ^^i'" XI / ^^^ ^(v,6)---^(T-n,,^fc)det(/C(r„,,6;r„^,,6/));^;,^i, 

fc=0 ■ ni = l-^ nk=l 

= det(l + /C^), (3.4) 

where Q{tj,^) = — X(s ,oo)(0 (J = 1)2,- ■■ ,m).(xj is the characteristic function.) The 
kernel /C is the extended Airy kernel, 

, ,7r <iAe-^'"-")AiK, + A)Aifc + A), n > r„ 



b) For f3o = (3- or (3o = (3+, we have JJ.^I ) with a different kernel. It is denoted as /C12 and 
is given by 

^12(7-1, 6; T"2, 6) 

/C2(ri, ei; r,, 6) - Ai(6) /„°° rfAe-^^Ai(6 + A) + Ai(6)e^-«-^ rs < 0. 

cj In the region where (3q < /?_ or (3q > /?+, the fluctuation is equivalent to those of the 
Brownian motion. In terms of the Fredholm representation, when (3i < [32 < ■ ■ ■ < Pm < 
/5_, we have 

lim F[Hi^~\f3u^.) < si, ■ • • ,hI^-\(3^,^_) < s^] = det(l + ICQ), (3.7) 

where the kernel is 



c? 



}Cg-{Pi,^uP2,^2)={ 



^''''^ A > P., 



v/27r(/3_-/3i)' 

2{/3_-/3i) 2(/32-/3i) 



«i te-ei)^ 



(3.^ 



^ y/27T(f3--l3i) ^27r(/32-/3i)' 



A</3: 



2- 



r/ie results for the case where P+ < Pi < P2 < ■ ■ ■ < Pm is analogous, 
(zz) When/3. > /3+. 

a)For (3 < f3c, the same fluctuation as ()3.8|) is obtained. 
(Hi) When (3- = (3+. 

a)For (3 < (3- and (3 > (3-, the same fluctuation as ()3.8|1 is obtained, 
b) For (3 = (3c, the results can be obtained as a limiting case of Theorem 5.1 in section\^ 

Remarks. 

1. The process characterized by the Fredholm determinant with the extended Airy ker- 
nel [30,31] in (i-a) is called the Airy process [21,22]. This is the same as the process 
of the largest eigenvalue in the Dyson's Brownian Motion model [29] between unitary 
ensembles. 

2. For (i-c) we can easily calculate the multi-point joint distributions. For instance the 



two-point joint distribution is calculated as 



lim P[/Jif-)(A,7_) < Si,/fif-)(/?2,7- 



= 1 



1 

+ 2 



d^iJCcAPu^uPu^i 



-) < ^2] 
d^2JCG.{P2,^2;P2,^2) 



+ 



00 ^oo 

1 J S2 

oo 



oo 



Si 



d^i / rf6 



S2 



Si 



d^i 



/CG.(/3i,6;/3i,ei) /CG_(/?i,ei;/?2,6) 

/CG.(/32,6;/3i,ei) /CG_(/?2,e2;/?2,e2) 

/Cg_(/32,6;/32,6) /CG_(/?2,6;A,ei) 

/CG„(/9i,6;/92,e2) /Cg_(/9i,6;/9i,6) 

ii2-iir 



,2(/3_-/32) 



S2 g 2(/32~l3i) 

OO ^V2^(/52-/5i)V2vr(/3--/32)' 



(3.9) 



Note that the 3 by 3 determinant, det[KG_(ti, ^i;tj, ^j)]fj^i becomes zero when tj = 
/?!, /32 since the 2 by 2 determinant at each Pj, det[KG_ {Pk, 6j Pk, ^j)]lj=i (^ = 1)2), 
vanishes. Hence bigger determinants also vanish. The integrand in ()3.9|) represents 
the propagation of a Brownian particle. 

Proof. Here we prove (i-a), (i-c), (ii-a) and (iii-a). The remainings, (i-b) and (iii-b), will 
be obtained as corollaries of the results in the following sections. 

First we start our proof with the fact that equal time multi-point correlation of height 
fluctuations at the odd time M = 2N — 1 has the Fredholm representation [22] , 

P [h{n, M) < h, h{r2, M),<k--- , h{r^, M) < /^] 

oo -. m m 

= 5Z ^ 5Z Xl ■ ■ ■ 5Z XI 9{ri,,xi) ■ ■ -giri^, Xk) det {KN{ri^,Xi] r^^,, XiO)f;z'=i 

fc=0 Ji=l xi ik=^ ^k 

= dei{l + KNg), (3.10) 



where 



girj,Xi) = -Xii,,oo){xi), 

KN{ri = 2ui,xi;r2 = 2^2, X2) = kN{2ui,xi;2u2,X2) - (p2ui,2u2{xi,X2), (3.11) 

Kn{2ui,xi;2u2,X2), 

(l-a)2(«2-«i) r dzi f dz2zl^ zi (1 - a/zi)^-i+"i(l - 0^2)^"^" 



-«2 



X 



(27rz)2 Jc^^ ^1 7c«, ^2 zTzi-Z2{l-az^)^-'-^-{l-a/z2Y-^+'-- 
1 - 7-/2:1 1 - 7+2:2 



1 - 7-/2:2 1 - 7+2:1 ' 



(l-a)^("2-''l) 



(t)ri,r2{Xl,X2) = { 2^* 



!ci ^Z''^-'=^[{1 - az){l - «/^)]«2-«i, U2 > Ml, 



U2 ^ Ml. 



(3.12) 



(3.13) 



In ()3.12|) . C/j. means a contour with a radius Ri and Ci is the unit circle. In both 
cases, they enclose the origin anticlockwise. One takes the radiuses Cr^, Cr,^ in a way that 
7_ < Cr^ < Cr^ < I/7+. Note that ()3.10|1 is valid for 7+7. < 1. 

Next we discuss the asymptotics by applying the saddle point method to the ker- 
nel (fTT^ and dSIini). 

proof of (i-a) 

We prove ()3.4j) and ()3.5j) along the same line as the derivation of the Proposition 4.1 in [25]. 
We give only the outline of the proof. 

• Asymptotics of K^ 

We set 

KNin = 2ui, xi, r2 = 2u2, X2) 

[^-(^)' 'I ^^1 / (^^2Z2 ,^(9p,.^i(.l)+9,2.^2(lA2)) 



{2my Jc,^ z, Jc,^ Z2 ^--^(/^-i) 

Zi- Z2I- 1-1 Z2 1 - 7+2:1 

where /3i = mi/A^, /?2 = —U2IN, /ii,/i2 are arbitrary constants at this stage and 

g^^p{z) = (1 + /3) log(^ -a) -{I -(3) log(l - az) -{^l + (3) \ogz. (3.15) 

Using a(/5), c(/3), (i(/3) defined in ()2.(jp . ()2.13|1 and ()2.14p respectively, we consider the scal- 
ing, 

A = /3o + ^r,, (3.16) 

Xi = a{Pi)N + d{l3o)NH,. (3.17) 

When we fix /Uc(/?) = a(/3) + 1, two critical points of g^^piz) are combined to the double 
critical point pdP) given by 



PciP) = p(/ic(/3), /3) = ^^^=^^^^S, (3.18) 

Vl - p + ay/1 + p 

where ^;^(^),^(Pc(/3)) = <,(/?), /3(Pc(/?)) = 0. Since 7_ < Pc(/5) < :^ for a < 7 < ^, 
P- < P < P+, we can deform the contour of Zi to 
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Then we find 



+ X2iPo)cm'N'^'r^ + XsiPoHPofr! + '^wl (3.21) 



Z2 



A2(/3o)c(/3o)'iVV3^| _ A3(/?o)c(/5o)V| + '-wl (3.22) 



where 



^i(/^o) = -^,9Mf^)APc{P)) 



(3.23) 



/3=/9o 



Similarly one gets 

Z\ - Z2 T2-T1+ l{Wi + W2) 

l,+Nii-,.iM) ~ (Pe(/?o))(«^-^^)'^(''")^'''e«--«-^+^«^-^+^«-^ (3.25) 

Substituting (prT!I|) - (IT^ into ^HM . one obtains 

gAi(/3o)c{/3o)Af2/3(ri-T2)+A2(/3o)c(/3o)27Vl/3(r2-r|)+A3(A))c{/3o)3(r3-r|)+6T2-eiTi 



^ ^ rfwi /" rfwa f -, r ) e'«i"'i+*«^"'^+t('"?+'"2) 



4vr2 Ji^^,=r,. """' Jlmwo=r,o ^^' \ T2 - n + i(wi + W2 



(3.26) 






e 

"OO 



rfAe-^("i-"2)Ai(ei + A)Ai(6 + A), (3.27) 

where r^j > is a convergence factor. 
• Asymptotics of (j)2ui,2u2{xi,X2) 

Next we consider the asymptotics of 02mi,2m2(^i?^2)- For ui < U2 

(Pri,T2\^li ^2} 

^ ' ' z"2-"4(l-az)(l-a/z)]"2-"i 



27ri ./ci z 



'^ "' ' _^^2-A^(Mc(/32)-l)-a;i+iV(/.c(/3i)-l)gAf9^,(/5i),,3iW+A^9Ma^2),/32(V2^). (3.28) 



27ri ./c^ -z 
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We set 



z = P.iPo) { 1 + ^^^) - vm (l - ^(^(-1 - -)) (3.29) 

^ ^l- ...!.n/. (^2-^a)V (3.30) 



Applying these to g^,,{|3^),|3^{z) and fl'^,(/32),/32(l/2;) respectively, we get 

^^Mc{/3i),/3i(^) ~ ^^A'cCA),/?! (Pc(A)) - -(ri - ior)^ (3.31 

Ng^lc{P2),|32 ( 7 ) ~ ^^Mc(/32),/32(Pc(/52)) + ^(^2 - i<y) 



(3.32) 



One also obtains 

^x2-7V(Me(/32)-i)-:r'i+7V(Mc{/3i)-i) _ (p^(/3Q))'^('^'')^''''(^2"^^)e'''(^2-Ci)_ (^3^33) 

Hence one finds 

v{6-€i)d(/3o)Af'/^ 



^3 ^3 



gAi(/3o)c(/3o)A^2/3(^^_^2)+^2{/3o)c{/3o)2Af^/='(r2-r|)+A3(/3o)c(/3o)3{r3-r|)-^ + ^ 

— / rfae^«^-«^+"i~"')"-(^2-^i)'^' (3.34) 

gAi(/3o)c(/3o)A'2/3(ri-r2)+A2{/3o)c(/3o)2A'^/^(r2-r|)+A3(/3o)c(A,)^(r3-T|)-:i + ^ 

/oo 
^^g-A(n-r2)Ai(^^ + A)Ai(e2 + A). (3.35) 

-oo 

From ()3.27|) and ()3.35p we get the desired expression. Note that the pre-factor of the 
extended Airy kernel does not affect the Fredholm determinant. 

proof of (i-c),(ii-a) and (iii-a) 

One can prove (i-c) in a fairly similar manner to (i-a) but an essential modification is 
necessary for the proof of the other two cases. This results from the fact that the model 
defined by ()2.ip - ()2.5|) is not well-defined in the latter case since the parameter of geometric 
distribution at i = 1, j = 1, which is equal to 7+7-, is greater than 1. Thus we consider a 
slightly modified model in which the parameter of the random variable at the point (1,1) 
is zero as in [27]. When 7+7- < 1, the modification is unnecessary but we consider the 
modified model here because the modification does not change the asymptotic properties 
of the model and allows us to treat all cases in a parallel fashion. 
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For 7+7- < 1 we can get easily the relation of the equal time multi-point correlation 
function between these models by generalizing the relation in one-point case [27]. Let 
h'^{ri,t) represent the height in the original model and h{ri,t) the height of the modified 
model. One has 

F[h{n, t = M = 2N-l)<h,--- , Krm, t = M)< Q 
^ P[/i+(ri, M) < Zi, ■ ■ ■ , /i+(r^, M) < U - 7+7-P[/i+(ri, M) < Zi - 1, ■ ■ ■ , /i+(r^, M) < l^ - 1] 

1 - 7+7_ 

(3.36) 

When 7+7- < 1, P[/i+(ri, M) < /i, ■ ■ ■ , h'^{rm, M) < lm\ is represented as the Fredholm 
determinant ()3.1()j) with the kernel 1)3.111) - ()3.13|1 . The problem is that the Fredholm 
determinant is not well-defined for 7+7_ > 1 (See (3.45) below). We would like to obtain 
another representation applicable for the case where 7+7- > 1 by modifying the Fredholm 
representation. 

We start from the Fredholm representation ()3.1()j) for 7+7- < 1. Using the relation 

Zi 1 - 7+2:2 1 - l-/zi ^2 , 1 - 7+7- 



2:1 - 2:2 1 - 7+^1 1 - 7-/2:2 Z1-Z2 (1 -7+2:1) (1 -7-/2:2)' 
the kernel Kn can be divided into two terms, 

kN{2Ui,Xi]2u2,X2) 

dzi f dz2 Z2^ Z2 F{ui,Zi) 



(3.37) 



{2niy Jc^^ zx Jc^^ Z2 zf zi - 2:2 F(m2, Z2) 



. 1 f dzi 1 ^, ,1 1 /" dz2 z'r, 



X2 



2m Jc^^ zi zl^ ' 1 - 7+2:1 2Tii Jc^^^ Z2 F{u2, Z2) 1 - 7-/2^2 

= K2Niui, xi] U2, X2) + (1 - 7+7_)£'+(mi, xi)D_{u2, X2), (3.38) 

where 



^(^' ^) = .. \,2un ' ,.w-i-. - (3.39) 



(1 - a/2)^-i+- 
;i-a)2«(l-«^) 

Remember that this expression was derived under the condition a < 7_ < i?2 < -Ri < 
1/7+ < 1/a. The Fredholm determinant can be deformed to 

det [1 + Knq] 



m 00 



= det[l-FK2Arfi'] < 1- (1 -7+7_) X^ X] X] D_{Uk,X2)E{Uk,Xi]Uk',X2)D+{Uk',X2) 
I k,k'=lxi=lk X2=lf.i 

(3.40) 
where 

K2n{.Ui,Xi]U2,X2) = K2n{Ui,Xi;U2,X2) - 4>2ui,2u2{Xi, X2) , 

E{ui,xi;u2,X2) = (1 - K2n)~^{ui,xi;u2,X2). (3.41) 
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In the above equation let us focus our attention to the contribution of terms, 

m oo 

(1 - 7+7_) Y, Yl ^-(^'^' ^)D+{uk, x), (3.42) 



k=l x=li. 



which arises from the 'delta function' part of E{ui,xi;u2,yi). Here one divides D± into 
two parts, 

1 f dz2 z% 1,7! 



L)_(m,x) = — - / — ^ + Tv z = Dx-^{u,x)^D2^{u,x), 

27rz Jc^, Z2 F{u, Z2) 1 - 7-/^2 F[u, 7_) 



R'r, 



D+Kx) = -i. / ^l-^^i!^ + 7^FKl/7+) = /^i+(«,x) + D2+Kx), (3.43) 
2vr2 Jc. zi zfl- 7+2i 

where the radiuses of contours R'^ are taken to satisfy a < i?2 < 7- < 1/7+ < -^1 < l/o^- 
Notice the second terms appear from the contribution of the poles at 7_ in D-E{u, x) and 
at 1/7+ in D^{k,x). One finds 



(1 - 7+7_) Y^ D+{uk, x)D_{uk, x) 

m 

-- (1 -7+7_)^^ \Di+{uk,x)Di^{uk,x) + Di+{uk,x)D2-{uk 



+ D2+{uk,x)Di_{uk,x) + D2+{uk,x)D2^{uk,x) . (3.44) 

The last term can be rewritten as 

(1-7+7-) E D2Au.x)D24u.x) = (1-7,7-) £(7+7-)^^^^^ = (7+7-)'^ ^^|^ 

(3.45) 
One notices that the series on the middle is divergent when 7+7- > 1 but that the 
difficulty is avoided in the right most expression. In addition one can easily find that 
such a difficulty does not arise for the other terms in (|3.44j) and the remaining terms 
of (1 — 'y+'y_)'^D_{uk,Xi)E{uk,Xi;Uk',X2)D^{uk',X2) in (j3.4Up with the contribution of 
fin^ subtracted, 

m 00 00 

E E E D^{uk,xi)E'{uk,xi;Uk',X2)D+{uk',X2), (3.46) 

k,k'=l xi=lk X2=ly 



where 



c,/. X ,E{uk,Xi]Uk',X2)-5l{xi,X2), for k = k, 

E{uk,Xi]Uk',X2) = < r , /,/ (3-47) 

E{uk,Xi;Uk',X2), tor fc 7^ fc , 
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and 

oo 

Y,St{x,y)f{y) = f{x). (3.48) 

y=ik 

Eventually one obtains the representation of F[h^{ri, M) < /i, ■ ■ ■ , h~^{rjn, M) < Im], 

{m cx) oo 

^ ^ ^ D^{Uk,Xi)E'{Uk,Xi]Uk',X2)D+{Uk',X2) 
k,k'=l xi=lk X2=ly 
m 

+ ^^ yDi+{uk,x)Di_{uk,x) + Di+{uk,x)D2-{uk,x) + D2+{uk,x)Di_{uk, 



(3.49) 



which is well-defined for 7+7- > 1 as well. Note that, when 7+7- > 1, the original meaning 
as a probability is lost in ()3.49p . 

Next we consider the asymptotics of ()3.49|) . We set Ui = PiNJi = aG_{Pi,'y-)N + 
dG{l-)N2Si 
ixi = aG_{f3i,'y-)N + dG{'y-)N2^A and A^ ^ 00. Let us notice 

D,4u,x)=^-^^^^ [ % \^ _e^....M^), (3.50) 

27r« Jc„, z [l-a z){l - -f+z) 

where u = N(3 and x = Ny. Since 

y = ac^ iP, -/-) + dci-f-)^ = yo + 6y, (3.51) 

N2 

the saddle point p^ for gy+i^/^{z) is known to be 

Changing the path of 2; in a way that it crosses p^, 

z = Pc\'^ + r = 7- 1 + ^^ r , 3.53 

'I i/3.- (3)dGN'2 / \ {(3-- P)dGN-^ ' 
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we get 



-- 9y+i,(3{Pciy,P)) + l^g'y+iAPciy^P)) i^-p^iv^P)}' 



= 9yo+i,p{Pciyo,P)) 
In the second equality we used 



c^G^logT- 



Ni 



2{(3o -f3)N' 



(3.54) 



9y+l,/3{Pciy,P)) ~ 9yo+l,0iPciyo,P)) + 

= 9yo+iAp'^{yo,P))- 



'^9y+l,l3 



dy 



5y + 



1 d'^gy^ 



l,/3 



y=yo 



2! d'^y 

e 



Sy' 



y=yo 



2{(3o-(3)N' 



(3.55) 



and 

1 



2! 



9y+iAPciy^f3)) {^-Pciy,f3)} ~ 7rA+iAPciyo,f3)) {^-^^(y,/?)} 



w 



2!- 



2{Po-P)N- 
(3.56) 



Combining the above, one finally finds 



D,+ {u,x)r^V,+ {P,0 



F(m,7_) 



'2(/3_-/3) 



(1 - 7+7-)7^ y/2n{P_-P)dGMm 



(3.57) 



Similarly, one can show 

Di_(w,x)~Pi_(AO 



{2 dc(7_)-^ 

-|^ g 2(/3-/3+) dG(7+)^ 

F(«, i/7+)(i - i+i-)Y+ y2H5^:^gg±|rfG(7-)ivi 



aG-(7-)-aG+(7+) „ dg(7^.) ^ 1 (ag (7- )-aG(T+)) ''G(T+) „Ar 

Xe (/3-/3+) ''dG(7-)'^" g 4(/3-/3+) dG(7-)'''\ (3.58) 



From these relations the asymptotics of K2N is also obtained easily, 

K2NiuuXi;u2,X2) ~ (1 -7+7_)Pl+(/5l,^l)Pl_(/32,6)• 



(3.59) 



In fact one can find that ()3.59|) does not contribute to the Fredholm determinant by the 
following discussion. First we consider the asymptotics when /5+ = /5_ = /3o deforming the 
contour of Zi such that they cross the double critical point pdPo), 



K2n{PoN,xi;PoN,X2] 



dzi 



dZ2zl'-''^^^-^^ 



(2vr02 Jc^^ z, Jc,^ Z2 ^--^(/^-i) 



,^(Smc,/3o(^i)+9mc,/3o(1/^2)). 



Z2 



Z\ - Z2 

(3.60) 
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where 



zi ~ Pc{I3q) 1 



IWi 



From (IT^ - (IT^ and noticing 



1 

Z2 



IW2 



PciPo) V c?(/3o)iVV3 



(3.61; 



z 



X2+N{l~ficil3o)) 



xi+N{l-fj.c(.M) 



PciPo 



,dGAf^(6-6) 



exp 



. ,dc 1 ac- — a 2 . .dc^ji^ , og- — a 2, 



(3.62) 



we obtain 



i^2^~(Pc(/?0))"^^'(«^-«^^ 



1 



rfA^l/3 

ciAAi(%iVlei + ^^^nI + X)M{^Nk, + ^^^-^ivl + A). (3.63) 



d 



d 



Using the asymptotics of Airy function 



d d 



1 e 3" d 



Ai(^^A^6^i + ^^- m + x) 



(3.64) 



where ac- — a > 0, we can find K2N goes to asymptotically. Combining this fact and 
the relation 

k2N{PoN,Xr,PoN,X2) 

~(l-7+7_)Pi+(/?o,6)^i-(/3o,6) 

1 (l-a/7_)i+*(l-«/7+)i-* 



;i - 7+7_) 



2 -I N 



(^^^_)«G-(7-) (1 _ a7_)i-/3o (1 - Q,^_^) 



l+/3o 



X e 4(/3o-/3+) dG(7_) 



one finds 






_(7+7„)"G-(7-) (1 _ a7_)i-^o (1 - a7+)i+^o 
Next we consider the case where Pi ^ ^2- From (|3.59|) . it is straightforward to see 

k2N{PlN,Xup2N,X2) 

~(l-7+7-)^i+(/?i,6)^i-(/32,6) 

~ (1 - 7+7-) X 7^^-^M(l - "/7-)(l - "7-)}"'"* 



(3.65) 



< 1. (3.66) 



X 



1 (1 - a/7_)i+^2 (1 - af-f+y-^^ _ K5(2-)zf^G(^£G(2+)" ^ 



(7+7_)"G-(7-) (1 - a7_)i-/32 (1 _ a7_^)i+. 



/32 



(3.67) 
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Hence from ()3.66|) and ()3.67p one finds that the contribution of K2n{PiN,xi; P2N,X2) is 
of order 0{e~^) and is neghgible in the Fredholm determinant of ()3.49|) . 

To obtain the asymptotics of ()3.49|) . one also needs the asymptotics of (f). Let us 
represent (f)ri,r2{xi,X2) as 

0.,,,(xi,x,) = ^^ ~ ^/!"""^ / ^e^/--(^), (3.68) 



27ri Jc^ z 

where 

fyuviiz) = 9yj+iMz) - 9y2+iM^)- (3-69) 

We scale yi as 

y, = aG-(A,7-) + ^%^ ^yo + 5m, (3.70) 

and adjust the path of z such that it crosses the saddle point p{. of fyi,y2iz), 

^ = Pc[^ + -, T^T^] ~ 7- I 1 + j-r^i^i -^2 + tw)]. (3.71) 



Then we get 



f (.\ f (.\ ^Glog7-(6-6) (6 - 6)' + W^ .^7.s 



From ^^^ and (jT7^ . one finds 



2 






dcN^ 



[I -a) 



(l-a7_)(l-a/7_) 



^iP^-P-) ewlV 



v/27r(/32-/3i)' 



(3.73) 



Substituting these asymptotic forms to ()3.49p and picking up the terms which do not 
vanish asymptotically, one finds 



P[/i+(ri,M)</i,---,/i+(r„,M)</^] 

F{uk, 1/7 



det[l + K2M9] { E(7+7-)'- ^^7'^^^;\ - 1 + A(/5, s)) 



fc=i 



m 



,^^4 v/27r(/?_-A) 



^ /-oo /-oo e 2('32-/3i) e 2{/3i-ft-i) e ^c/s.-ft) 

l^ik.Z:'<J^H 4, v/2vr(/?2 - A) V27r(A-A-i)v/2vr(/5_-A) 

(3.74) 
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where (7-|-7_)''' ^"^' '^\' K{(3, s) represents the contribution of the terms included in ^ D^E'D^ 
such as 

k,k' x\,X2 

y^ y^ -D2-(wfc,a;i)02«fe,2«fc.(a;i,a;2)-ftr2Ar(ufc",a;3;Ufc',a;2)-D2+(Mfc',a;2), 

k,k" ,k' xi,X2,x-j, 

y^ y^ ^2-(^fc,a;i)02tifc,2«fe..(a;i,a;3)</>2«fe..,2«fc/(a;3,a;2)-D2+(^',a;2), (3.75) 

k,k" ,k' xi,X2,X3 

and so on. In fact the first term in ()3.74j) cancels due to the subtraction in (|3.36j) : 

det[l + ir27v^] 



det[l + i^2^,] 1^(7^7-)'^ ^|^^i^( - 1 + A(/?, .)) I 

U=i J^yuk,!-) J 



7+7- det[l + K^m-i] \ Efa^-)'^"' ^F?'^^^!^ ( - 1 + ^(^' ^)) 



-^TTT E (t^ det [1 + K2.g]) ^^7' ^^^j^ ( - 1 + A(A s)) (7.7^)^ 
y_)A^2 ^ V^Sfc / ^(«fc,7-) 



~ — - — ; — r Y". ^2N{uk, X = ac-N + dcN^Sk] Uk, x = aa-N + dcN^Sk)—^ — \- 

(-l + A(/3,5))(7+7-)''= 
~ 0, (3.76) 

where gi_i{rj,Xj) = ~Xii ■-i,oo){x)- Thus when we consider the asymptotic behavior 
of (jSSni), P[/i+(ri, M), ■ ■ ■ ',h+{rJ,M)] can be replaced with 



r[/i+(ri,M),---,/i+(r„,M)] 

e 

^2n{P^-Pk) 



1 - E / ^^^ 



fe=l -"^fe 



-{ei-?i_i)2 



00 ^00 



^ h k,h<k. 4i ' ' ' 4, "^^^ ' " '^^V2vr(/?2 - Ay ' ' ' v/2vr(A-A-i)V2vr(/3--A)' 

(3.77) 
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Using the fact that ()3.36|) can be expressed as 

P[/i(ri,M)</i,---,/i(r„,M)</J 

1 ™ o \ 

~ r [/i+(ri, M) < /i, ■ ■ ■ , /i+(r^, M) < /^], (3.78) 

one finally gets 



P[/i(ri,M)</i,--- , /i(r^, M) < Z„] 

1 - E / ^^'^ 



fc=l ■'^fe 



e 

"t — , 

' v/2vr(/?_-A) 



-(g»-Si-l)^ 



i^ik.t"<J^H 4. V27r(/?2 - A) V27r(A-A-i)v/2vr(/3--A) 

(3.79) 

This expression is the same as the Fredholm representation in ()3.7|) with the kernel ()3.8p . 

D 

We end this chapter by providing a proof of (ii-a) in Theorem 2.1, for which the same 
strategy as the above proof of (i-c),(ii-a) and (iii-a) is applicable. The only difference is 
the asymptotics of Di_ ()3.58|1 . This changes to 

g2 de(7_)2 

I g 2{/3c-/3+) dG(7+)2 



Di_(M, = /3,iV,x)~P'i_(/5c,e) 



F(«„ l/7+)(l - 7+7-)7^ V27r(/3-/3+)^da(7-)iV^ 



da (7- ) 

(3.80) 



This means that we can obtain the Gaussian as a scaling limit for both Di+(n,x) and 
Di-{u,x) since the (3c point is the crossing point of two lines with Gaussian fiuctuation. 
Together with this, the asymptotics of K2N also change; 

K2n{uc, xi; Uc, X2) ~ (1 - 7+7-)^i+(/5c, ^iWUPc, 6) ~ Oie-""). (3.81) 

Considering these two relations in addition to ()3.57j) and ()3.73p . we pick up the terms 
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which do not vanish in ()3.49p . 

P[/i+(rc = 2Me,M) </] 



det[l + gK,^] { {^^^^)^l^^plM. (_i + A(;3^, ,)) 

+ 1-/ ^6 . .^ =- / ^6- 



V2vr(/5- - Pc) Js V27r(/3e-/3+)gg4 



^0(7-) 



JI 'ic(7_)^ 



+ / c^^i , ^ / d^2 , . , , , (3.82) 

Js V2vr(/3- - /5.) y. V2vr(/5c-/5+)g{^ 

Here we can neglect the first term since this term vanishes due to the subtraction in ()3.36p . 
Hence one obtains, 

.„ , f^ e 2(/3--/3c) fS p 2(/3c-/3+) dc{7+)2 

^— J-00 V27r(/?_ - /?,) J-00 v^27r(/?,-/?+)gn±2 

This completes the proof of (ii-a) in Theorem 2.1. In principle, one can also apply the 
same method to the case of the multipoint function including the (3c point. 

4 Transition around GOE^ 

4.1 Limiting Kernel 

Let us suppose that, when 7_ = 70, the limit shapes of a(/3) and ac-iP, 7-) cross at /3 = Pq. 
We call this point the GOE^ point. 

In this section, we obtain the kernel describing the multi-point height fluctuation near 
this GOE^ point. We consider the case where the parameter a < 7+ < I/70 is fixed and 

0/ GrtQldG nVci 



7_ scales like 



^-=^°('-rf(^^)' ^'-'^ 



with u fixed. The result is 
Theorem 4.1. 

lim P[i/7v(ri, Po) < si, ■ ■ ■ , HN{Tm, Po) < s.^] = det(l + KG), (4.2) 

TV^oo 

where g{Tj, ^ = -X(s,,oo)(0 6' = !> 2, ■ ■ ■ , m) and 
^(n,^i;T2,6) 
'/C2(ri,a;T-2,6)+Ai(6)/o"rfAe-(-+-^)^Ai(6-A), ^ ^ uj + T2>0, 

./C2(ri, ei; T2, 6) - Ai(ei) /o°° rfAe(-+-^)^Ai(6 + A) + Ai(ei)e%^-«^ (-+-), u; + r^ < 0. 

(4.3) 
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Remark. This kernel seems to be new. If we set u; = from the beginning, this gives 
Theorem 3.1 (i-b). Notice that when we focus on the one-point correlation function, the 
case where r = 0, u; = r' is essentially the same as the case where t = t',uj = 0. 

Proof. First we give the contour integral representation of the kernel. For 7_ < i?2, it 
reads 

k^iri = 2ui,xi;r2 = 2u2,X2) 
(1 - a)2(«2-"i) f dzi f dz2Z^'' zi (l-a/2i)^-i+"i(l-a22)^~^""' 



(27r^)2 Jc^^ zi Jc^^ Z2 zl^ z^- Z2{l-aziY-^--^{l-a/z2)''-^+^- 

1 -7-/^1 1 -1+Z2 
1 -7_/z2l -1+Zi 

For 7_ > _R2, we have to add the contribution of the pole at 2:2 = 7_, resulting in 

k^iri = 2ui,xi;r2 = 2u2,X2) 
(l-a)2(«2-«i) r ^^^ r ^^^ zf Zi (1 - a/zi)^-i+"i(l - a^2)^"^~"' 



X 



(4.4) 



(27rz)2 Jc^^ zi Jc^^ Z2 zl'zi-Z2{l-aziY-^-^^{l-a/z2Y-'+^- 
1 - 7^/2:1 1 - 7+Z2 
1 - 7_/2;2 1 - 7+Zi 

(1 - a)2(«2-"i) f dzi 1 (1 - a/2i)^-i+"i7!2 (l-a7_)^-i-"2(^l-7+7_) 
Cn, ^Wil-az.Y-^-'^'il-l+z,) (1 - a/7_)^-i+"2 ' 

(4.5) 



Now we consider the asymptotics when we set (jHH), 

X, = iVa(/5o) + rf(/3o)iV'/'e., (4.6) 

r^ = 2u, = 2N (/5o + ^) , (4.7) 

with 2 = 1,2 and take N ^ 00. The analysis is almost the same as in [25] and hence the 
details are omitted. When u + T2 > 0, one uses ()4.4p and deforms the contours of zi and 
Z2 such that 

^i-7o(l + ^^^^(n-.^i)), (4.8) 

---(^- -T7W^Wmi^2 + tW2)) . (4.9) 
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Then we get 



^ ^Xi(l3o)c{l3o)N'^/HT^-T2)+\2{Mc{f^ofN^/^T^~Ti)+X3{l3o)c{l3of{Tf-Ti)+^2T2-iiiTi 

X — n / dwi / (iW2 I r H 

4vr^ Jlm«,i=r,i Jlmw2=V2 \ r2 - Ti + l{Wi + W2) UJ + T2 + IW2 

X g«6«'i+«6«'2+§(w5'+u.|) (4.10) 

with the definition of Aj given in (j3.23|) . Noticing 

— - / dwi / (iw2 I H 

47r^ yim«,i=r,i Jlmw2=V2 V ^3 - Ti + ?(Wi + W2) CU + Ts + 2^2 

= / e'^("i-"^)Ai(ei + A)Ai(e2 + A)(iA + Ai(ei) / rfAe-("+^2)^Ai(e2 - A), (4.11) 
Jo io 

one gets the desired expression for this case. 

The UJ + T2 < case can also be treated in a similar fashion. In this case we consider 
the asymptotics using ()4.5p since the contour of Z2 ()4.9|) cross the real axis on the left of 
7o. The second term in the RHS of ()4.5p produces the third term in the RHS of ()4.3|) . D 

4.2 Connection to the Baik-Rains analysis 

When we specialize to the case of the one-point height fiuctuation, our formula reduces to 

lim ¥[Hn{0, /3o) <s] = det(l + JCQ), (4.12) 

where Q{^) = -X(s,oo)(0 and 

/C(x, y) = IC2{x, y) + A{x)B{y, cu), (4.13) 

with 



K,2{x,y)= Ai(x + A)Ai(y + A)t/A, (4.14) 
Jo 

A{x) = Ai(x), (4.15) 

{/ rfAe-'^^Ai(x-A), u > (4.16a) 

- / dXe^^Aiix + A) + e^'/3-x..^ ^ ^ q_ (A.Wh) 

Notice that the expression on the right hand side of ()4.12|1 is independent of /?o and that, 
for the special case where /3o = 0, another expression for the same quantity was previously 
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obtained in [27]. Hence the two expressions should be the same. For uj = 0, this was 
already shown in [32]. In this subsection, we prove the equivalence for any u generalizing 
the arguments in [32]. 
We first proceed as 

det(l + ICg) = det(l + (/C2 + A ® B)g) 

= det(l + IC2Q) det(l + (1 + /C2^)-M ® BQ) 

00 POD 



= F2{s)\l-J dxj dyp{x,y)A{x)B{y,uj)], (4.17) 

where p{x, y) is the kernel of the operator (1 + /C2^)^^- Hence if one defines 

dx I dyp{x,y)A{x)B{y,uj), (4-18) 

J s 

one has 



Let us also define 



and 



det(l + ICg) = F2is)a{s, iu) = F{s). (4.19) 

b{s,uj)= / dypis,y)Biy,uj), (4.20) 



Q{x)= / dyp{x,y)A{y), (4.21) 

J s 

q{s)=Q{s). (4.22) 

Now we show that the functions a, b have the following properties. 
Proposition 4.2. 

—a = qb, (4.23) 

d 

7-6 = qa- ujb, (4.24) 

OS 

d 

-_a = q^a - (g' + uq)b, (4.25) 

OUJ 

d 

7— & = (q' - ujq)a + (cu^ - s - g^)6. (4.26) 

OUJ 
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Proof. From the results in [5], one has 

m \- f \ x+f \ Q{x)P{y) - P{x)Q{y) 

R{x,y) =p{x,y)-6^{x-y) = , 4.27 

x-y 

§-^Q{y) = q{s){S^{y -s)- p{s, y)), (4.28) 

Q\y) = P{y) - u{s)Q{y) + q{s)p{s, y) - q{s)5+{y - s), (4.30) 

P'{y) = yQ{y) - 2Q{yHs) + u{s)P{y) + p{s)p{s, y) - v{s)b^{y - s), (4.3i) 

u^ -1v = g^ (4.32) 

4 = V — Qu, (4.33) 



where 



p{x,y)Ai'iy)dy, (4.34) 

p{s)=P{s), (4.35) 



u{s) = / Q{x)Ai{x)dx, (4.36) 

v{s) = / P{x)Ai{x)dx, (4.37) 



S+{y-s)f{y)dy = f{s). (4.38) 

J s 

Besides, we also use 

^B{y, uj) = Ai{y) - ujB{y, u), (4.39) 

■^B{y, uj) = Ai'iy) - ujAi{y) + {uj' - y)B{y, uj). (4.40) 

These equations can be shown immediately. For the case where u; > we use ()4.1(ia|l and 
compute 

—Biy,u) = -j d\e-'XAi{y-\) 

POO 

= -/ dXe--^{yAiiy-X)-Ai'\y-X)}. (4.42) 

Jo 

In ()4.42|1 the Airy equation, Ai"(x) = xAi(x), is used. These equations lead to ()4.39|1 
and ()4.40|1 respectively. For the case where uj < they can be shown by applying the same 
method to (I4.16bll . 
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Now the first two equalities, ()4.23|) and ()4.24|) . can be shown as 

d f°° d 

—a = Q{s)B{s,uj)- dy—Q{y)B{y,(jj) 

/■oo 

= qis) dyp{s,y)B{y,uj) = qb, (4.43) 



d f^ f d \ f°° 

h = -j dy[—p{s,y)\B{y,uj)-q{s) dyQ{y)B{y,u) - p{s,s)B{s,uj) 

/•oo g 

= qa-q+ / dyp{s, y) g-B{y, u) 

= qa- cub. (4.44) 

For the third equality ()4.25|1 . one starts from 

d f°° d 

—a = - Q{y)—B{y,uj)dy 

poo poo 

= -/ Qiy)Ai'iy)dy + uj Q{y)M{y)dy 

J s J s 

/•CXD /»00 

-cuM Q{y)B{y,Lu)dy+ yQ{y)B{y,uj)dy. (4.45) 

J s J s 

The second and the third terms are easily seen to be uju and — ti;^(l — a) respectively. 
Using ()4.30|) . the first term is computed as 



"OO 

1 oo 



Q{y)A^{y)dy = - [Q{y)Ai{y)] ^ + / Q\y)Ai{y)dy 

J S 

/OO /'OO /"OO 

P{y)Ai{y)dy-u Q{y)Ai{y)dy + q p{s,y)Ai{y)dy 
J s J s 

= v-u'^ + q^. (4.46) 

Next we consider the last term. First we calculate 

yQ{y) = P\y) - uQ\y) - Q{y){u' - 2v) + (p - 5+){qu - p) 

= P\y) - uQ\y) - Q{y)q^ - pg' + 5+g', (4.47) 

where we use fl4.3Up . ()4.31|) in the first equality and ()4.32|) . ()4.33|) in the second equality. 
Then one can show that the last term is 

/»00 /"OO /"OO 

/ dyyQ{y)B{y, uj) = -q^ / dyQiy)Biy, w) -q'b+ dyP\y)B{y) (4.48) 

O S v S rJ S 

/oo 
dyQ\y)B{y, w) + q'B{s, w) 

= q^a - {q' + ujq)b - v + u^ - q^ - uju + uj'^il - a), (4.49) 

26 



where we use fj4.39|) and ()4.40|) . Combining these, one gets ()4.25|) . 
Finally, for the proof of the last equality ()4.26|) . one starts from 

d f°° d 

dyp{s,y)ki{y) 

/•OO /"OO 

+ uj^ I dyp{s,y)B{s,uj)- dyyp{s,y)B{s,y). (4.50) 

J s J s 

The first three terms are easily seen to be p = g' + uq, —ujq, uj% respectively. The last term 
can also be calculated. The term yp in the integrand can be shown as 

l/p(-5, y) = {.P- qu)Q{y) - qQ'iv) + (g^ + s)p{s, y) - q'^6+{s - y) (4.51) 

due to (H?f|l and (IQUl) . Thus from (IOH|) . (IOn|) and (H3T| . one gets 

P-OO P-OQ 

- dyyp{s,y)B{s,y) = -{q -uq) dyQ{y)B{y) - {q^ + s)h - uq (4.52) 

J s J s 

These equations lead to ()4.26p . D 

This proposition, combined with a(s, 0) = 6(s, 0) = e~-'s '?(^)'^^ proved in [32], shows 
that the functions a, 6 as defined by ()4.18|) . ()4.20j) are the same as those in [27] with the 
identification a-Q^ = a, 6br = —b, u = —q, w = 2uj. Thus the one-point height fluctuation 
defined in (I4.19p corresponds to that in [27] under the above identification. 

4.3 GOE^ to GUE/Gaussian Transition 

Using the above results, it is possible to study the fluctuation properties of the PNG model 
quite in detail. In this subsection, we set c<j = and consider the transition from GOE^ 
(r = 0) to GUE (r -^ oo) and Gaussian (r -^ — oo). As an example, we here consider the 
average of the scaled height. For the case where r = or r = oo, we can easily compute 
the average numerically using the Painleve expression of the GOE^ or GUE [5,6]. For the 
other values of r, however, such a representation has not been known. To compute the 
value numerically, one can use the differential equations of a and 6, (J4.23|) - (j4.26p . On the 
other hand, we also obtain the asymptotic behaviors of the average for each case where 
r ~ 0, oo, — oo as follows. 

• r ~ 

The equations ()4.23|) - ()4.26p enable us to know the behaviors of a for r ~ 0. We get 



a(s, r) = a(s, 0) + 



da{s, t) 



dr 



T + 



1 d^a{s,0) 



2 <9r2 

T = U 



27 



2 1 d^a{s,T) 



T^ + 



3! dr^ 

T=U 



r^ + . . . . (4.53) 

T=0 



hen the distribution behaves as Substituting ()4.53|) to F{s) defined in ()4.19|) . the average 
behaves as 

/ sF'{s)ds ~ / sF^Q^2{s)ds + t / ei{s)ds + ^ e2{s)ds + ^ / e^{s)ds -\ 

= -0.49364 - 0.89941r + 0.41582r^ - 0.12409r=^ + ■ ■ ■ . (4.54) 

Here ej(s)'s are computed as 

e,{s) = se-^^'^^^^'^ {P2{s){q{sf - q{s)') - F,{s) {q{sf + sq{s) - q\s)q{s))] , (4.55) 

e^is) = se-/^°°''(^)'^^ [P^is] {q{s) (-g(s) + q{sY + 2q{sfs + s^ -sq'{s) - g'(s)')) 

+F2(s) {-q{s) + q{sf + s q\s) - g'(s)') } , (4.56) 

e3(s) = se-^»°°^(^')'^" {P2{s) {q{sf + q{sf + 2q{s) s - {q{s) + q{sf + q{sf s + s^) q\s) 

+ {-qisr + s)q'isf + q'isf) 

+ F2{s) {q{s) (2 - q{sf - q{sf + q{s) s - 3 g(s)^ s - 3 q{sf s^ - 5^+ 

(qis) + qisr + qisY s + s') q'{s) + {qisf + 2 .) q'{sf - q'isf) ) } , 

(4.57) 

using (jOni)- (P3K|1 and a{s, 0) = b{s, 0) = e-i'"'?(^)^^ proved in [32]. 

• r — i> oo 

One can also get the asymptotic behavior of a as r — > oo. By noting 



1 /"°° 
B{x, t) = - dee-^A\{x - 0/t) 



r r^ r 

one finds 

r>00 /"OO 



M£)_^ + ^_.... (4.58) 



1 /'°° /"°° ■u(s) 
a(x,r)~l / dx j dy p{x, y) Ai{x) Ai{y) = 1 ^—^, (4.59) 

where u{s) defined in (j4.36p reads 

qixfdx. (4.60) 

Then one has 

/ sF'{s)ds ~ / sF^{s)ds - - / F2(s)u(s)cis 1.77109 + -. (4.61) 
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• r ^ — cxD 

Finally we consider the asymptotics where r ~ — oo. Due to ()4.23|) . the probability density 
is represented as 

^ = q{s)b{s, t)F,{s) + a{s, t)F^{s). (4.62) 

We evaluate the first term of the above equation as follows. From the results of Baik- Rains, 

lim a{s, t) = 1, 

T— >00 

lim b{s, r) = 0, 
a{s,T) = b{s,-T)e^'-'% (4.63) 



one finds when r ^ oo, 

Let us scale the variable s as 



^3 



b{s,T)^e^ '^. (AM) 



s = T^y, (4.65) 

and we assume y > 0. (This can be justified in the below discussion.) Using the asymptotic 
behavior of the Airy function, we also find 

g(.) ~ Ai(.) ~ ^-^e-i^l (4.66) 

Then due to ()4.64p . ()4.66p and F2{s) ~ 1 as s — i> oo, one finds 

q{s)b{s, r)F,is) ~ ^^y^^^TJ^^^ '''\ (4-67) 

where 

9iy) = ly'^' -y + i- (4.68) 

The condition g'{y) = leads to y = 1. Rescaling y as y = 1 + y'/kP''^ ^^^ expanding 
g{y) around y = 1, one gets 

qis)b{s, T)F,is) ~ ^-^i^e^. (4.69) 

On the other hand, for the second term, one finds 

a{s, t)F^{s) = a{r^y, T)F^{T^y) ~ (4.70) 

since a{s,T) ~ 1,^2(5) ^ as s -^ 00. Thus the probability density is approximated by 
the Gaussian with the peak at s = r^ when — r is large sufficiently and one gets 

f sF2is)ds r^ t\ (4.71) 
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Using ()4.54|) , ()4.6H) and ()4.7H) . we can calculate the average of the distribution numer- 
ically. This is shown in Fig 3. In principle, other statistical quantities such as two point 
function analyzed in [33] can be also studied. 

Before closing the section, we would like to mention that the GOE^ distribution seems 
rather universal [34] but the numerical values of statistical quantities of the GOE^ have 
not been given in the literature. Since we have the data at hand, we list them here. The 
known data for GOE/GUE/GSE [36] are also presented for comparison. (In this table, 
the data for GSE is shown according to the notation in [18]. Therefore the values of the 
average and standard deviation in GSE are different from those values written in [35]. ) 





average 


s.d. 


skewness 


kurtosis 


G0E2 


-0.49364 


1.1100 


0.3917 


0.309 


GOE 


-1.20653 


1.2680 


0.2935 


0.165 


GUE 


-1.77109 


0.9018 


0.2241 


0.093 


GSE 


-3.26243 


1.0176 


0.1655 


0.049 



5 Transition around Fn 



5.1 Limiting Kernel 

Let us suppose that, when 7_ = 70 = I/7+, the limit shape is tangent to that in bulk (J2J 
at P = jSq. Let us call this point the Fq point. In this case we again have to consider the 
modified model as in the proof of (i-c), (ii-a) and (iii-a) in Theorem 3.1 since the parameter 
of geometric distribution at the {i = l,j = 1) point is close to unity in the original model. 
The relation between the modified and original models is given in ()3.36p . 



P[/i(ri, t = 2N)<h,--- , hir^, t = 2N) < I 
-- P[/i+(ri, 2N) < /i, ■ ■ ■ , /i+(r^, 2N) < /J + 

-P[/i+(ri,2A^) </i-l,- 



7+7- 



1 - 7+7_ 
/i+(r^,2iV)</^-l]}. 



{P[/i+(ri,2iV)</i,---,/i+(r^,2iV)</, 

(5.1) 



When we consider the case where 7+7. = 1 we always treat P[/i+(ri, t) < /i, ■ ■ ■ , h^ir^, t) < 
lm\ in the right hand side as the Fredholm determinant ()3.10|) with the kernel ()3.11j) - ()3.13|) . 



In the scaling limit where 



Ti = 2(3oN + 2riA^2/3 



li = a{ po + 



Ari/3 



N + d{Po)N-^s„ 



7± 



1 

7o 



UJ± 



d{(3o)N^/' 



(5.2) 



30 



one finds 

lim P[ifAr(ri,/3o) < Si, ■ ■ ■ , ifAr(r„,/3o) < s. 



N^oo 



1 ™ a \ 

1 + ^ Et^ J™ P[i^^(ri,/5o)<Si,---,iJ+(r„,/?o)<sJ. (5.3) 






We obtain the kernel describing the multi-point height fluctuation near this Fq point. 
The result is 

Theorem 5.1. 

When CJ+ + u;_ > 0, 

(1 "^ r) \ 

1 + ■ V ■— det(l + }cg), 
00^ + iO^ ^ OSj j 

(5.4) 
where Q{tj,C) = -X(s,,oo)(0 6' = 1> 2, ■ ■ ■ , m), 

^(n, Ci; ^2, 6) = ^2(ri, 6; ^2, 6) + (^+ + cu-)S(c^+, n, ^i)B'{w_, ra, 6), (5-5) 

r/-rfAe-(-+--)^Ai(6-A), c.+ -ri>0, 

i3(u;+,ri,4i) = < -r3+^i (5.6) 

\-/,°°rfAe(-+--^)^Ai(ei + A) + e^^-«^(-+--^), ^+ - n < 0, 

r r rfAe-("-+-2)^Ai(e2 - A), u^ + T2>0, 



/o°° rfAe('^-+"2)^Ai(e2 + A) + e^i ^(^-+^2)^ uj_ + T2<0 



Remarks. 



1. The Fredholm representation ()5.4j) and ()5.5j] includes the terms which are not analytic 
for u+ + uj- < 0, i.e., the difficulty of divergence that we discussed in the proof of 
Theorem 3.1 (i-c), (ii-a) and (iii-a) arises. However, one can avoid this in a similar 
manner as for the Theorem 3.1. 

For example let us consider the term included in this Fredholm determinant which 
corresponds to the last term in ()3.44|) . 

POO 



(cu+ + c^_)e"++"- / rfee"^^"++"-^ (5.8) 

J s 

For u;+ + u;_ > 0, this can be calculated as 

31 



Here though the LHS is meaningful only for c<j+ + tu- > 0, the RHS is analytic for 
all region of ojj^ + uj_. Hence the difficulty of divergence of (jS.Sj) can be avoided if 
one replaces the integral by the RHS in ()5.9j) when uj+ +uj- < 0. Other terms which 
diverge when u;+ + u;_ < can be treated in the same manner. 

2. When we set uj± = we obtain Theorem 3-1 (iii). 

Proof. For the case where uj+ — ti > 0,uj_ + T2 > 0, the contour integral representation 
of the kernel reads 

Kn[2Ui,Xi]2U2,X2) = - 



X ^-^^^-^-^-l^y (5.10) 

1 - 7+2:1 1 - 7_/Z2 

The asymptotics can be studied in the same manner as in the previous section with the 
result, 

^ gAi(/3o)c(/3o)Ar2/3(ri-r2)+A2(/3o)c(/3o)''Af'''^{Ti2-r|)+A3(/3o)c(/3o)='{Tf-r|)+6r2-6n 

X -^ /" dwi I dw2e'^''"'+'^^'"^+-M+-^l) 

1 UJ+ + U^ 



T2 - Ti + i{Wi + W2) (CU+ - n + iWi) {OJ- + T2 + ^^2) 

'(5.11) 

The one arrives at the desired expression. 

In the case of a;+ — ri > 0,a;_ + r2 < the contribution of the pole at Z2 = 7- are 
added as 

(l-a)2(«2-«i) r dzi 1 (1 - «/;zi)^-i+"i7!^ (l-a7_)^-i-"2(l-7+7„) 
^i Jcn, ~^ (1 - ^^^'^^'-"^(l - 7+^1) (1 - a/7-)^-i+"2 

^ (1 _ «)2(«2-«i)(p^(^^))(6-5i)</3o)iVV3_ 



^ gAi(/3o)c(/3o)Af2/3(ri-r2)+A2(A))c(A))2Afl/3(r2-T|)+A3(/3o)c{/3o)3(rf-r|)+C2r2-Cin 



X e^ 6(r2+<^-)^ / dwi—^ (5.12) 

27r Jlrnw^=m CJ+ - Ti + IWi 

(|5.11|1 and (J5.12J) lead to the desired expression. Similar to this case we add the term of 
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pole at zi = 1/7+ in the case of u^ — ti < 0, u^ + T2 > 0. This term reads 

(1 _ a)2(«2-«i) r dz2 ,, (1 - az2)'^"'~"^7? (1 - a7+)^-i+"i(l - 7+7„) 



27ri Jc^^ Z2 ' (l-a/z2)^-i+"^ (7+^2-1) (1 - a/7+)^-i-"i 

U aj [p,[Po)) rf(/5o)iVV3 

-r3+c^3 , . 1 /■ p«6«'2 + |«'i 

X e^r^-«i(-+--i)— / dw2— (5.13) 

27r Jlmw2=V2 UJ^+T2 + IW2 

Finally we consider the case of u;+ — ri < 0, cj- + r2 < 0. In addition to (JS.III) . ()5.12|) 
and (j5.13p we have to consider the term caused by the product of the residues, 

(^ — n"v \N-l+ui /I _ \N-1-U2 

^ ^ ^+ ^~ (1 - a/7+)^-i-"i (1 - a/7_)^-i+"2 

y. gAi(/3o)c{/3o)Af2/3(ri-r2)+A2(/3o)c(/3o)2Afl''^{T2-r|)+A3{/3o)c{/3o)^(rf-r|)+6r2-Cin 



'^^"^ -6(T2+t^-)^ \ + -6(^+-n) 



X e^ ?2(r2+c._jg 5 52(0.+ -nj^ ^5_^4^) 

Combining each term we get the desired expression. D 

5.2 Connection to the Baik-Rains analysis 

When we speciahze to the case of the one-point height fluctuation, our formula again should 
be equivalent to the ones obtained in [27]. We show the equivalence in this subsection. 
When we set Tj = 0, the Fredholm determinant reads 

dei{l + ICg) = F2{s) h - {uj+ + uj_) j dx dyp{x,y)B{x,uj+)B{y,uj_)\, (5.15) 

where F2{s) is the GUE Tracy- Widom distribution and B{x,u) is defined in ()4.16p . We 
have the following proposition for this expression. 

Proposition 5.2. 



CO poo 



l — {uj^ + uj-) I dx I dyp{x,y)B{x,uj+)B{y,uj-) = a{s,uj+)a{s,uj-) — b{s,uj+)b{s,uj-). 

(5.16) 
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Remark. Using the proposition and ()5.3|) . we can easily calculate the one-point height 
fluctuation. 

lim F[Hn{t = 0, Po) < s] 

N—*oo 

1 d \ 

IH ■ — {a{s,u+)a{s,iJ^) - b{s,uj+)b{s,uj_)]F2{s) 

tU-j- + UJ— OS J 

I 

= a(s,iu^)a(s,uj_)F2(s) H \a(s,uj^)a(s,u) — b(s,u^)b(s,uj)}F2(s). (5.17) 

0;+ + u;_ 

This corresponds to the expression of the one-point fluctuation in [27]. 

Proof. Let / denote the function on the left hand side. By differentiation, one finds 

^/(s) = {u+ + u_){a{s,uj+)a{s,u.) - f{s)}. (5.18) 

This can be read as 



oo ^oo 



s 

oo 



f{s) = a{s,u+)a{s,u_) + — I dxl dyp{x,y)B{x,uj+)B{y,ujJ). (5.19) 

The second term is now computed as 

— dx dyp{x,y)B{x,uj+)B{y,uj_) 

/OO /'OO O /"OO 

dyp{s,y)B{y,uj^) + dxB{x,uj+)— dyp{x,y)B{y,uj_) 

POO 

= -B{s,uj+) dyp{s,y)B{y,uj-) 

J S 
/•OO /"OO 

+ / dxB{x,uj+) dy{- p{x,s) +6+{x- s)}p{y,s)B{y,u_) 

/oo 
dyp{y,s)B{y,u^) 

= -b{s,uj+)b{s,U-). (5.20) 

where we use the equation (2.16) in [5]. D 

6 Conclusion 

We have studied the PNG model with external sources. We have represented the multi- 
point equal time correlation functions of the height fluctuation as the Fredholm determi- 
nant. In order to get these quantities we have first identified the limit shape of the model. 
The results are summarized in Theorem 2.1. The limit shape consists of a straight line and 
a circular one. There are three cases according to the value of 7+7_ and they are illustrated 
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in Fig. 2. The statistics of the one-point height fluctuation is Gaussian on the straight hne 
and the GUE Tracy- Widom distribution on the circular line. In addition, there are special 
points denoted by /3± where the two limit shapes meet and their statistics of the height 
fluctuation is GOE^ when 7+7- < 1 or Fq when 7+7- = 1. 

Our main results are summarized as Theorems 3.1, 4.1 and 5.1. The spatial configura- 
tion of the height fluctuation over some region of the surface has been shown to converge to 
the stochastic processes represented by the Fredholm determinants obtained in the above 
theorems. In addition to the one-dimensional Brownian motion and the Airy process 
which describe the correlations near edges and in the bulk respectively, we have obtained 
the Fredholm determinant which describes the correlation near the points (3±. 

These quantities vary with the parameter of the nucleation at edges 7± in fj2.4|) and ()2.5p . 
For the case of 7+7- < 1 the Fredholm determinant obtained in Theorem 4.1 describes 
the transition between GOE^ and GUE/Gaussian. On the other hand, when 7+7^ = 1 the 
points I3± merge to /3o and we have obtained the Fredholm determinant expression with 
the transition kernel between Fq and Gaussian. These Fredholm determinants have the 
parameters Tj which determine the position in the surface. When we set r^ = r 7^ we can 
obtain the one-point distribution at the point away from (3±. Fig. 4 shows the comparison 
of these transitive distributions with the Monte-Carlo simulations of the PNG model in 
the corresponding points. These figures confirm our analyses. 
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Figure Captions 



Fig.l: Rules of the discrete PNG model. Fig. (a) shows an example of the rules 1 and 2. 
A step with height k at the origin is generated according to the geometric distribution. It 
grows laterally one step in both directions in the next time step. 

Fig. (b) shows a collision of the steps with height one and two. The height of the 
origin, which is the point they collide, is two by the rule 3. 



Fig. 2: Typical examples of the limit shape of the PNG model with external sources. The 
solid lines show the snapshots of Monte Carlo simulations for A^ = 5000, a = 0.32. The 
dashed lines represent the limit shape in bulk given by (j2.6j) . The parameters 7+ and 7_ 
are 7+ = 0.79, 7_ = 0.63 in (a), 7+ = 1.58, 7_ = 1.26 in (b) and 7+ = 1.58, 7_ = 0.63 in 
(c). 

Fig. (a) corresponds to the region (i) in Theorem 2.1 where the hight fluctuation in 
bulk obeys the GUE Tracy- Widom distribution while the fluctuations near both edges are 
described by Gaussian. Fig.(b) corresponds to the region (ii) where the height fluctuations 
at all points obey the Gaussian distribution. Fig.(c) represent the case of region (iii). The 
limit shape is a straight line tangent to the limit shape in bulk. 



Fig. 3: Average in the GOE^ to GUE/Gaussian transition. Three lines represent the 
asymptotics in the case where r ~ 00, 0, — cxd respectively. + gives the value obtained by 
solving (|4.23|) - (|4.26|) numerically. 



Fig4: The limiting distribution near the competing points between bulk and edge. In 
Fig. (a) the three curves are the distributions for r = 1,0,-1 from the left obtained in 
Theorem 4.1. Each point is the data of a Monte-Carlo simulation near the origin where 
t = 2000(A^ = 1000), a = 0.32,7+ = 0.32, 7_ = 1.0 and x = 352,0,-352 respectively. 
Fig.(b) shows the distribution near Fq. The curves represent the transitive distributions 
for r = 0, 1 from the left using Theorem5.1. The points are numerical data for t = 
2000(A^ = 1000), a = 0.32,7+ = 1.0, 7_ = 1.0 and x = 0,352 respectively. 
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